clear; clc;
close all

% 均匀半空间电场衰减情况
mu = 4e-7 * pi;
S = 0.1;
fre = 10;
% fre = 1;

Omega = 2 * pi * fre;
k = sqrt(-sqrt(-1) * Omega * mu * S);
z = 0:10:10000;
Ex = exp(-k * z);
plot(real(Ex), -z/1000, 'r');
hold on
plot(imag(Ex), -z/1000, '--');
xlabel('E_x/E_0');
ylabel('深度/km');
legend('实部值', '虚部值');